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Abstract 

Color transparency CT depends on the formation of a wavepacket of small spatial 
extent. It is useful to interpret experimental searches for CT with a multiple scattering 
scattering series based on wavepacket-nucleon scattering instead of the standard one 
using nucleon-nucleon scattering. We develop several new techniques which are valid 
for differing ranges of energy. These techniques are applied to verify some early approx- 
imations; study new forms of the wave-packet-nucleon interaction; examine effects of 
treating wave packets of non-zero size; and predict the production of N*'s in electron 
scattering experiments. 
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I. INTRODUCTION 

Color Transparency (CT) is the postulated [1,2] absence of final (or initial) state 
interactions caused by the cancellation of color fields of a system of quarks and gluons 
with small spatial separation. For example, suppose an electron impinges on a nucleus 
knocking out a proton at high momentum transfer. The consequence of color trans- 
parency is that there is no exponential loss of flux as the ejected particle propagates 
through the nucleus. Thus, the usually "black" nucleus becomes transparent. We re- 
strict our attention to processes for which the fundamental reaction is elastic, or at 
least a two-body reaction. This requires that the nuclear excitation energy be known 
well enough to ensure that no extra pions are created. This subject is under active 
experimental investigation [3-5]. 

The existence of color transparency depends on: (1) the formation of a small-sized 
wave-packet in a high momentum transfer reaction. (2) the interaction between such a 
small object and nucleons being suppressed (color neutrality or screening) and (3) the 
wave-packet escaping the nucleus while still small [1,2]. That color neutrality (screening) 
causes the cross section of small-sized color singlet configurations with hadrons to be 
small was found in Ref. [6], and is well-reviewed in Refs. [7-9]. So we take item (2) 
as given. The truth of item (1), for experimentally available energies, is an interesting 
issue. It is discussed in Refs. [9-11], and is not probed in depth here. 

It is also true that at experimentally available energies, the small object does expand 
as it moves through the nucleus. Thus the final state interactions are suppressed but 
not zero. The importance of this expansion was found by Farrar et al. [12], and by 
Jennings and Miller [13]. See also Ref. [14]. The purpose of this paper is to determine 
improved methods to calculate the evolution of a wave packet as it moves through a 
nucleus. In Ref. [12] the uncertainty principle is used to argue that the wave-packet- 
nucleon interaction can be treated by using an effective cross section, cr e ff(Z) oc Z 
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where Z is the propagation length. In Ref. [13], the wave-packet-nucleon interactions 
are computed using a hadronic basis. An approximation to the full multiple scattering 
series is made in which the first-order term is exponentiated. The work of Ref. [14] 
uses an exact Green's function for the evolution of the wave-packet in a model in which 
the wave-packet-nucleon interaction is proportional to the sum of the squares of the 
transverse separations between the quarks in the ejected wave-packet, b 2 . Moreover, 
the baryon states are treated as quarks bound in a two-dimensional harmonic oscillator. 
The exact solution depends on using this particular interaction and model space. 

The aim of the present work is to develop techniques and approximations that allow 
the use of more general models of the wavepacket-nucleon interaction and baryon space. 
A general formal expression for the necessary multiple scattering series is derived in 
Section II. Different approximation schemes are also defined in that section. Numerical 
studies of the different approximations are made, using the (e, e'p) reaction as an exam- 
ple, in Sect III. (The initial wave packet is assumed to be of zero size.) Approximations 
that work best in different regions of energy are identified. The applications in Sect. Ill 
result from using a wave-packet-nucleon interaction proportional to b 2 . Sect. IV shows 
how different interactions, e.g. that of Ref. [15], can be used. This demonstrates the 
general nature of our methods. Sect. V displays how our techniques are used in situa- 
tions in which the wavepacket is allowed to have a non-zero size. Another application 
is to N* production, Sect. VI. The final section is reserved for a few brief summary 
remarks. 



II. FORMALISM AND APPROXIMATIONS 

Our technique is to treat the wavepacket, formed in high momentum transfer pro- 
cesses, as a coherent sum of baryon states. The means that the relevant wave equation 
should be a matrix equation. Hence the starting point for our analysis is the relativistic 
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wave equation 

(V 2 +p 2 )\V N , p ) = U\y N , p ). (2.1) 

The quantity \^n, p ) represents a vector in internal quark and nuclear (center of mass 
position) spaces. The subscripts refer to the asymptotic boundary condition that a free 
nucleon (N) has a momentum p. The operator p is 

p= \J p 2 + M 2 — M 2 , (2.2) 

M 2 -M 2 , s 

« P+ ^— , (2.3) 

where the energy eigenstates of the quark "Hamiltonian" , M 2 , are such that 

M 2 \n) = \n}M 2 . (2.4) 

In this notation, the nucleon is the ground state, with n = N with eigenvalue M 2 . The 
matrix operator U acts in both quark and nuclear space. It has matrix elements between 
the nucleon and its excited states and also depends on the nuclear density p{R). 

We wish to simplify Eq. (2.1). To this end, let R be the operator denoting the 
position of the center of mass of the struck nucleon. Also, we take the final nucleon 
momentum, p, to be large and in the Z direction. Define another column vector *& p (R) 
as 

* p {it) = (it\* Njt ). (2.5) 

It is convenient to factor out a plane wave factor and write, 

* p (R) = e i * z <p(R)\N). (2.6) 

This is the defining equation for <p(R), a matrix in baryon space. Then, as an operator 
equation, assuming the momentum is large compared with any gradient in the problem, 

2zp^^ = U(R)^(R) : (2.7) 



where 

U (R) = e~ l P z U{R) e^ z . (2.8) 
The solution to Eq. (2.7) is a path ordered exponential, 

<f(R) = V exp [j ^U(B, . (2.9) 

Here, V is the path ordering symbol. Notice that we have chosen to use outgoing 
boundary conditions and that R = B+ZZ. We often suppress the index B in evaluating 
the path integral. 

Before proceeding it is worthwhile to discuss the nuclear interaction operator U (R) . 
This is the product of the wave-packet-nucleon interaction 4nf by the target nuclear 
density p(R) 

U(R) = 4Trfp(R). (2.10) 

In the forward scattering approximation used here, the operator / acts on the internal 
quark space only. We write / = f(b 2 ) where fry is the transverse (perpendicular to p) 
separation operator between quarks i and j in the wave packet and b 2 = b 2 j . That the 

ij 

interaction does not depend on the longitudinal coordinates is suggested by the work 
of Ref. [6]; see also Ref. [8]. Although there is no detailed theory for /, some general 
things are known. For small wave packets with b 2 <ti b 2 H where b 2 H = (N\b 2 \N), f should 
vanish. Interactions do occur for larger wavepackets. The expectation value of / in a 
normal baryon should be the experimentally measured forward scattering amplitude. 
Such are dominated by their imaginary part so we keep only Imf. Thus we write 

4tt/ = -z2^/(6 2 ), (2.11a) 

where a is the proton-nucleon total cross section with 

fib 2 = 0) = 0, (2.116) 
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(N\f(b 2 )\N) = l. (2.11c) 
In the two gluon exchange approximation, 

Bm m » J. (2.1U0 

We shall use a general function f{b 2 ), subject only to the constraints of Eqs. (2.11). 
The presence of the factor p allows the high energy baryon-nuclear cross section to be 
independent of the beam energy. Replacing p by p in Eq. (2.8) is of the same order of 
approximation as those already made and produces no essential changes in our results. 

To be specific, we examine the (e, e'p) reaction [4] and compute amplitudes for 
processes in which a proton is knocked out of some shell model bound state, a. With 
semi-exclusive kinematics, as in the experiment done at SLAC, we look for an exiting 
proton which has momentum, p, very close to that of the virtual photon, whose momen- 
tum we label by q. By close we mean that no extra pions are produced in the collision. 
Then, we define the "knockout" amplitude and cross section by the overlap 

M a = (N,a\T H (Q 2 )\V N , p ), (2.12a) 

aoc^|M a | 2 , (2.126) 

a 

where the sum is over all the occupied shells a. The cross section of Eq. (2.12) is 
obtained by integrating over the angle of the outgoing proton, corresponding to the 
SLAC experiment. We shall be concerned with ratios a/a B in which a B is obtained by 
replacing (R\^n jP ) of Eq. (2.12), as given by Eqs. (2.5), (2.6), and (2.9), by a proton 
plane wave function. 

The operator T H (Q 2 ) accounts for the absorption of a high momentum photon 
by a nucleon. Ignoring the effects of spin, the matrix elements of the hard scattering 
operator are 

(R, N\ T H (Q 2 )\R, m) = e~^ z 6& (R - R')F N , m (Q 2 ), (2.13) 
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where F NyTn (Q 2 ) is the inelastic transition form factor. The term e tqZ FN,m(Q 2 ) {m\ 
represents a high momentum wave packet propagating through the nucleus. 

After inserting a complete set of nuclear position states, the scattering amplitude 
M a has the expression: 

M a = j d s R$ a (R) (N\T H (Q 2 )e^ z V exp ^ ^ dZ 1 V(B,Z 1 ) ) j \N), (2.14) 

where 3> a (R) is the wavefunction of a nucleon bound to the nucleus in shell ct, and q 
is the three momentum of the incoming photon. We take the shell model parameter 
Tluj = 41 MeV. The operator V is given by 

V(R) = ^U(R), (2.15a) 

V(R) = -e~^ z - p(R) f{b 2 ) e^ z . (2.156) 

The "hat" over V is meant to symbolize that V(R) also depends on b 2 . We note that 
in general, V(R), V(R') ^ since the operators p and b 2 do not commute. All of 
the nuclear interactions are contained in the scattering wave function \^N,p)- At initial 
stages the wave packet is small and so are interactions. If the wave packet expands, 
the influence of interactions is more important. Thus, an accurate approximation must 
apply for both small and large sized wave packets. The aim of this work is to develop ap- 
proximation techniques for Eq. (2.9) or Eq. (2.14). This involves handling the evolution 
of an initially small wave packet as it moves through the nucleus. 

At this stage there are several things that one can do to obtain (p and We will 
consider what we call the "Order By Order" calculation, the "Low Energy Expansion" , 
then the "High Energy Expansion" and the "Equal Spaced Momentum" replacement, 
and finally the "Exponential Approximation" . 



A. The Order by Order Calculation (OBO) 
This means that we simply expand the exponential in Eq. (2.9) order by order. 

Vexpff dZ 1 V(Z 1 )\ = 1 + I dZ 1 V(Z 1 )+ I dZ 1 V(Z 1 ) f 1 dZ 2 V(Z 2 ) 

v-oo J — oo J — oo J — oo 

+ / dZ 1 V(Z 1 ) [ 1 dZ 2 V(Z 2 ) [ 2 dZ 3 V(Z 3 ) 

J — oo J — co J — oo 

+ ■■■ . (2.16) 

This approximation scheme should be accurate when V acts as a small number. This 
occurs at high enough energies such that wave packet expansion effects are minimal. 

B. The Low Energy Expansion (LEE) 
Before considering the low energy expansion in detail, we first obtain a general 
result. We write V as a sum of two terms, 

V(B, Z) = V L (B, Z) + AV L (B, Z), (2.17) 

with 

V L = (N\V\N), (2.18a) 

AV L = V-V L . (2.186) 

Notice that with the above definitions, (N\AVl\N) = 0. Further, if Vl and V are local 
operators in nuclear space, as in Eq. (2.15b), then 

[V L ,AV L } = 0. (2.19) 

Let us suppress the transverse nuclear coordinate, B, since it is the same for all terms. 
Then, in general, 



( Z A7 Sz dZ lVL ( Zl ) j Z * d Zl V L ( Zl ) 

+ / dZ 2 e Jz z AVL{Z 2 )e J - 00 

J — oo 



+ I dzj 1 dzj* dxVL{x) AV L ( Zl )e^ dxVL(x) AV L (Z 2 )ef-~ ^ 

J — oo J — oo 

+ ■■■ . (2.20) 
The low energy expansion is obtained by using Eq. (2.19) to simplify the above equation: 



p/.-.^W,) = W 1+ /* JZ^VUZO 

\ J — oo 



+ 1 dZ 1 AV L (Z 1 ) dZ 2 AV L (Z 2 ) + ■ ■ ■ 1(2.21) 



where 

X(B,Z)= / dZ 1 V L (B,Z 1 ). (2.22) 

J — oo 

This is called a low energy expansion, since at low energies the wave packet size increases 
and V « Vl- Then, AVl is small and one recovers the results of low energy nuclear 
physics. 

C. The High Energy Expansion (HEE) 
Here we also define a large Vh and small AVh potentials and treat AVh with 
perturbation theory. The separation of terms is now based on AVh acting as small at 
high energies. Thus we define 

V = V H +AV H , (2.23a) 
AV H = e^ z V e~^ z - V. (2.236) 

This approximation is expected to be valid at high energies because, in the closure limit, 
all momenta are equal and AVh = 0. The operators Vh and AVh do not commute, 
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so we must use the more difficult form, Eq. (2.20), of the path ordered exponential in 
Eqs. (2.9) and (2.14). 

D. The Exponential Approximation (EA) 
This is an approximation first made by Miller and Jennings which turns out to be 

useful to examine. These authors make the following approximation: 

-z 



\n) (n\e^ z V exp / V {Z 1 )dZ 1 \N) 

J — OO 



e ipZ exp 



\ n )( n \ / dZ 1 V{Z 1 ) 



\N). (2.24) 



At high energies V is small and the exponential is well approximated by its first-order 
expansion. The approximation is also accurate at low energies because the exponential 
accounts for an exact treatment of the usual distorted wave function. 

III. MODEL EVALUATIONS 
The aim here is to determine the regions of accuracy for the various approximations 
of Sect. II. We do this by using tractable models for the baryon space interaction 
operator V and initial wave packet. The baryonic states are chosen as those of a two 
dimensional (transverse) harmonic oscillator that binds two quarks. We take %uj = 
500 MeV . The interaction is specified as 

f(b 2 ) = b 2 /b 2 H , (3.1) 

in the notation of Eqs. (2.11). This is a simple function with the correct general prop- 
erties. A more elaborate form is studied in Sect. IV. Now consider the initial wave 
packet. We assume that the hard scattering, represented by T H , creates a wavepacket 
of zero-size. Such an object is called a Point-Like configuration, or PLC. Writing down 
only the internal (quark) degrees of freedom, we have 

(N\T H (Q 2 ) = F Njn (Q 2 )(n\, (3.2a) 

n 

= C(Q 2 )(b = 0\. (3.26) 
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So, F N:U (Q 2 ) = (b = 0\n)C(Q 2 ). For the two dimensional oscillator, the matrix element 
(b = 0\n) has the same value for all n, ^^^7^)' anc ^ a ^ f° rm factors are equal in this 
model. 

We now turn to the evaluation of Eq. (2.14). Inserting a complete set of oscillator 
states gives 

M a = f d 3 Re-* qZ $ a (R) jr F N!2m (Q 2 )e ip ^ z (2m\Ve W ( [ dZ 1 V(B,Z 1 )) \N) 

J m=0 \J-oo ) 

(3.3) 

where only the isotropic state n = 2m need be kept. The quantum number m represents 
the number of nodes in the 6-space baryon wavefunctions. With these preliminaries out 
of the way, we turn to the evaluation of the different approximations. 

A. Order by Order Calculation 
The order by order calculation is to evaluate Eq. (2.16). The result is 

M a = F NN J d 3 R^ a (R)e-^ z ^(R), (3.4) 

where 

^(R) = e^J2^ B °(R), (3-5) 

n=0 

and the index n is the order of the term in V. We find 

* OBO (E) = 1, (3.6) 



*? BO (fl) = -- j Z dZ lP (B, Z x ) (l - e ^-v){z- Zl )\ 5 (3 7) 

2 J-oo ^ ' 



^ BO (R) =(^Y [ dZ lP (B, Z 1 ) [ 1 dZ 2 p(B, Z 2 ) \l - e^-rt^O 

_|_ e i(p 2 -p)(Z 1 -Z 2 ) _ 2 e i(p 2 -p)(Z-Z 2 ) + 2 e i(P2-p)(Z 1 -Z 2 )+i(p 4 -p)(Z-Z 1 ) 



,(3.8) 
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*$ BO (R) = - (|) 3 / dZ 1 p(B,Z 1 ) [ 1 dZ 2 p(B,Z 2 ) [ 2 dZ 3 p(B,Z 3 ) 

J-oo J-oo J — oo 



X 



1 _ e i(P2-p)(Z-Z 1 ) + e i( P2 -p)(Z 1 -Z 2 ) _ 3 e i( P 2-p)(Z-Z 2 ) 

_|_ 2 e i (P2-p)(Zl-Z 2 )+i(p4-p)(Z-Zl) ^ ei(p2 _ p)(Z2 _ Z3) _ ei(p2 _ p )( Z _ Zl+Z2 _ Z3 ) 

+ Se^-P)^ 1 "^) - 9 e «(P2-p)(^-^3) + Q e i( P2 -p)(z 1 -z 3 )+i(p 4 ,-p)(z-z 1 ) 

_ 4 e i(p2-p)(^-^ 3 )+i(p4-p 2 )(^l-^2) _|_ 10e i(p4-p)(^-^ 2 )+i(p2-p)(^2-^3) 



Q e i(p6-p)(Z-Z 1 )+i(p 4 -p)(Z 1 -Z 2 )+i(p 2 -p)(Z 2 -Z 3 ) 



(3.9) 



where p is the momentum of the outgoing nucleon and p 2m = p 2 + M 2 — M 2m . Note that 
the above wavefunctions, ^^ BO , ^l/f? 50 , and ^l/^- 80 all vanish in the closure limit. The 
order by order curves are shown in Fig. 1. The curves labelled by OBO n are defined by 

n 

using q° B °™ = £ ^OBO 
m=0 

The distorted wave Born approximation DWBA is obtained by using a scattering 
wave function approximated by the neglect of excited baryon states. Then, since there 
are no quark space operators, the exact answer is obtained by simply exponentiating 
the first order result without excited states, * = e ipZ (l - § ff^ dZ 1 p(B, Z{^j. That 
is, 

Vdwba = e» z e-*f-- dZ ^^ = e» z ^ 3 &> (3.10) 

where 



a 



z 



\(B,Z) = -- dZ 1 p(B,Z 1 ). (3.11) 

^ J-oo 

The name DWBA arises from the notion that the plane wave factor of Eq. (3.10) is 
"distorted" by a factor, here an exponential damping factor. The DWBA curve is 
denoted by the straight solid line on all graphs. We see that for — Q 2 ~ 30GeV 2 the 
second and third order calculations are converging. At lower energies, there is poor 
convergence. This is expected, however, since at low energies we expect there to be 
many (more than two or three) scatterings. 
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B. The Low Energy Expansion (LEE) 
Now let us proceed to describe the low energy expansion. With the present b 2 scattering 
operator, the operators of Eqs. (2.17) and (2.18) become 

V L (B,Z) = ~p(B,Z), (3.12) 

AV L (B, Z) = ~p(B, Z)e~ ipZ Ab 2 e^ z . (3.13) 

where Ab 2 = 4 1. Here too, we have gone to third order in AVl- Using previous 

°H 

notation, we obtain 

^ EE {R) = e^ z \ (3.14) 



^ EE (R) = : [ Z dZ l p(B,Z 1 )e i ^- p ^ z - z ^, (3.15) 

2 J-oo 

z z 
^ EE (R) = (^Ye X{Z) [ dZ 1 p(B,Z 1 ) [ 1 dZ 2 p(B,Z 2 )\e i ^- p ^- z ^ 

^2/ J-oo J-oo 

— 2e i ( p2 ~ p ^ z ~ z ^ + 2e i ( p2 ~ p ^ Zl ~ Z2 ^ +i ( p4 ~ p ^ z ~ Zl ^ , (3.16) 



*% EE (R) = (^) 3 e A(Z) / dZ 1 p(B,Z 1 ) [ 1 dZ 2 p(B,Z 2 ) [ * dZ 3 p(B,Z 3 ) 

J-oo J-oo J-oo 

x Li(p2-p)(Z-Z 1 +Z 2 -Z 3 ) _ 2 e i(P2-p)(Z 1 -Z 3 ) + 4e i( P2 -p)(Z-Z 3 ) 

_ 4:e i(P4-p)(Z-Z 1 )+i( P2 - P )(Z 1 -Z 3 ) _ 8e i(p4-p)(Z-Z 2 )+ l ( P 2-p)(Z 2 -Z 3 ) 

_|_ 4 :e i(P2-p)(Z-Z 1 +Z 2 -Z 3 )+i( P4 - P )(Z 1 -Z2) 



+ Qe i ( p6 ~ p ^ z ~ Zl ^ +i ( p4 ~ p ^ Zl ~ Z2 ^ +i ( p2 ~ p ^ Z2 ~ Z3 ^ > 



(3.17) 



In the closure limit, the wavefunctions given by Eqs. (3.14 - 3.17) do not vanish but 
rather conspire to form the first three terms in the expansion of e~ x ( z \ times e +x ( z \ 
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Therefore, to all orders, the full wavefunction converges to the Born wavefunction, corre- 
sponding to transparency. Plots of cross sections using the above distorted waves appear 

n 

on Fig. 2. The curves labelled by LEE n are defined by using \& LEE " = ^ ^^ EE . 

m=0 

The DWBA curve is * = e ipZ ^ LEE . As before, the Born term is given by ^(R) = e ipZ . 
The fact that the LEE does not join smoothly up to the DWBA is evidence that the 
m = 1 "Roper" has influence below threshhold. 

C. The High Energy Expansion (HEE) 
In the present model, the evaluation of Eqs. (2.23a) and (2.23b) leads to 
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V H (B,Z) = ~p(B,Z)-^, (3.18) 



AV H (B,Z) = ~p(B, Z) 



- b 2 - b 2 



(3.19) 



Because Vh and AVh no longer commute, we must use the path ordered exponential 
shown in Eq. (2.20). The result is 

tff EE (R) = e* pZ fa i^fXf i{P2m ~ P)Z + J2 xn+1 ( Z y (P2m ~ P)Z ( 2m I £T I 2n ) 

_° [ Z dZl ^ Zl ) - Y e i{P^-p)Z e i{p 2n - P2l )Z, x n, z s 

li_o. 1 l-A(Z 1 )f- ; 1 v 



m,n,l 



x (2m | e ] 21) (21 \ ^ \ 2n) , (3.20) 



b 2 



where A(Z) is defined in Eq. (3.11) above, 



X(Z,Z 1 ) = ~ ( dZ 2 p(Z 2 ), (3.21) 
= ^fy, (3.22) 



and 



<*» i i a > ) § (")(!)? ) ■ < 323 > 
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The result (3.23) is obtained from 

Although it is not immediately obvious, ^^ EE (R), given by Eq. (3.20), in the closure 
limit, does indeed approach unity and color transparency is obtained. Notice that we 

— > 

have suppressed the transverse nuclear distance, B, in A, x, and \1/ for clarity. 

At this stage, a further approximation suggests itself. At very high energies, con- 
sidering the form of p in Eq. (2.3), it seems reasonable to replace 

P2m ~ P -> m(p 2 - p). (3.25) 

This is what we call the "Equal Spaced Momentum" (ESP) replacement. With the 
above replacement, we can now do the sum on / in Eq. (3.20) with the result that 

*™ - ^ ( + (*(*>- "Ml 1 -»(*>) 

1 ( ' l(l-A(Z))(l-y(Z,Z)) + ^. y{Z , Z) y 

2j_J Zl l-MZ 1 )^ n e X{Zl) 
((2n+l)(2m | e ^ | 2n> 

- (n + iyCP-P*)Si ( 2m | e MZ ' Zl) k \2n + 2) 

- ne l{p2 - p)Zl (2m \ e^^^ \ 2n - 2))^ . (3.26) 
The new function y is given by 

y(Z,Z 1 )=x(Z 1 )e^-^ z . (3.27) 

Intermediate states corresponding to m = n = 9 must be included to accurately obtain 
numerical results for expressions (3.20) and (3.26). On Fig. 3, we see this first order 
term of the HEE (solid) compared with the OBO s (dotdashed). As expected, the two 
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curves agree at higher energies. The ESP replacement is very accurate at all but the 
very lowest energies. 

D. The Exponential Approximation (EA) 

If one looks at only the first order term in Section IIIA, we see that the only thing 
which has changed between this wave function calculation and the traditional Glauber 
calculation is that the cross section a has been replaced by an effective cross section 
cr(Z, Z\) = cr(l — e t ( p2 ~ p ^ z ~ Zl ' ) ). Then, it is quite tempting to simply exponentiate this 
first order result. This so-called exponential approximation is not correct, since only the 
second excited state enters whereas we know that higher states enter in higher orders. 
However, because we have learned that the ESP replacement is a safe one, we can cast 
all higher excitations in terms of excitations of the second state. In fact, if we examine 
the difference A 2 between the exact second order OBO calculation and the second order 
term of the EA, we find that, 

z z 

A 2 = (°y[ dZ 1 p(B,Z 1 ) [ ^Z 2 p(B,Z 2 )e^-^- z ^(l-e 2 ^-^ z - z ^)\ 

J -oo J -oo ^ ' 

(3.28) 

where we have used the ESP replacement. Therefore, we see why the EA works as well 
as it does. The exponentiation of the first OBO contains many more of the higher order 
terms than one would naively think. Thus, at high energies, the difference between the 
EA and the true answer is rapidly suppressed, as shown above. The EA (dotted) is 
shown along with the first order OBO calculation on Fig. 4. 

E. Summary of Approximation Schemes 

We have presented many different approximations and graphs. Here, we compare 
and contrast the results of the above calculations. In particular, we want to see which 
approximation method or methods work best. To do this, we study the convergence 
properties among all curves. Looking at the Figures we see that: 
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1. The curves 0B0 2 and OBO3 converge at about — Q 2 30GeV 2 . This leads us to 
believe that above — Q 2 xs 30GeV 2 , the second/third order result is the true answer 
in this model. At low energy, we see no convergence. This is not surprising since 
at low energies, we expect there to be many (more than two or three) scatterings. 

2. The curves LEE 2 and LEE 3 follow each other closely up until -Q 2 ss 20GeV 2 . 
This indicates that LEE 3 is the correct answer at least until — Q 2 xs 20GeV 2 , and 
probably further, depending on how big the fourth order contribution is. Compar- 
ing Figs. 1 and 2, we see that LEE 3 and OB0 2 ,3 line up at about — Q 2 m 26GeV 2 . 
This gives us confidence that LEE3 is the correct answer for — Q 2 « 26 — 35GeV 2 . 
Therefore, we can conclude that LEE 3 is the correct answer all the way up to 
— Q 2 ~ 35GeV 2 . Of course, at higher energies, we expect LEE3 to be inaccurate, 
as did LEE\, and higher order terms would be needed. Not bad, however, for a 
low energy expansion. 

3. Fig. 3 shows the results of the high energy expansion including the ESP replace- 
ment. Also shown for comparison is the OBO3. The two curves are shown to line 
up for energies higher than — Q® 2QGeV® indicating that the true answer in 
this model is well approximated by the BEE. 

4. Fig. 4 shows the exponential approximation and the OBO\ for comparison. The 
EA matches the BEE all the way to -Q 2 « 6GeV 2 , indicating that the EA is a 
simple and fast way to calculate the distorted wave for high energies. 

5. The most successful of these approximations are the LEE 3 , the BEE and the EA. 
They are compared in Fig. 5. 

IV. MORE GENERAL INTERACTIONS 
The advantage of the present formalism is its applicability to general forms of the in- 
teraction, f(b 2 ) of Eqs. (2.11). The tests of Sect. Ill were performed with f(b 2 ) = b 2 /b 2 H . 
Here we use the form of Nikolaev and Zakharov [15], which is well approximated [16] 
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by 



fib") 



1 + 7 



1 _ e -^/b% 



(4.1) 



7 



This is derived for qq color singlet wave packets, but we use it here as a general repre- 
sentative of Eqs. (2.11). The parameter 7 = 0.762. 

We use the LEE to determine the cross sections predicted by this interaction. Then 
the result Eq. (2.18a) is unchanged and the new form of AVl is 



A straightforward application leads to the results of Fig. 6. At high energies the OB0 2 
and LEE 2 approximations are in agreement. At low energies the results of LEE 2 agree 
with those of the EA. Thus the LEE 2 provides a good representation of the correct 
answer for this model. We see that the results are qualitatively similar to those of 
Ref. [13], but the predicted ratios of <j/<jb are reduced by about 10%. 



So far we have assumed that a PLC was made in the hard scattering. Another 
consequence of this assumption, also used, is that the elastic and inelastic form factors, 
F NN and F 2rri)N , are equal. Yet another way to restate Eqs. (3.2) is to assert that 



nucleon. We now examine the effects of including a nonzero size for the PLC. In this 
section, we use the simpler form (~ b 2 ) for the wavepacket-nucleon interaction. 

We have incorporated the effects of finite size in two different ways. In the first 
case, we examine the effect which letting Tu(b) have a nonzero range has on the limit 
of — . Here, we assume 



AV L = ~p(R)- h-( 7 + l) e -^ 2 /^" . 

Z 7 L J 



(4.2) 



V. FINITE SIZE EFFECTS 



Tuib) oc 5^ (6), where b is the transverse distance between the two quarks in the struck 



T H (b) 



1 




(5.1) 



A 2 7T 
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It is then simple to calculate the form factors. In fact, a calculation similar to that for 
deriving Eq. (3.24) enables us to calculate all of the form factors at once. The result is 

F * m ' N = W^W) (a^t%) ' (5 ' 2) 

where we take bn = Ifm. In order to illustrate the effects of finite size, we examine 
the OBO\ case for simplicity. In the case of unequal form factors, the distorted OBO\ 
wave becomes 

V OB °i(R) = l- a -f dZ 1 p(B,Z 1 )(l- Ejb!L e *<J»-p)V-*)\ (5.3) 
2J-oo ^ Fn,n ' 

This is to be compared with Eqs. (3.6) and (3.7). Fig. 7 shows the ratio of a to ob for 
a very large value of — Q 2 . As A 2 approaches zero, Tn(b) approaches a delta function, 
and the transparency is unity. As A 2 goes to infinity, the ratio of the form factors goes 
to zero, and the first order term of the standard DWBA results. This shows that the 
predicted cross sections are controlled by the existence of a PLC. 

The immediate question is, then, what is A 2 ? Above, it is an arbitrary parameter. 
But, if we believe that a small sized configuration is produced in the hard interaction, 
then A 2 should be related to the momentum transfer of the incoming photon to the 
ejectile. Suppose we have a nucleon of three quarks which absorbs a photon of three 
momentum q. Each quark, then, gets approximately q/3 of the photon's momentum. 
So, we use a form [17] for the hard scattering operator 

T H (b)=e^. (5.4) 

This form is suggested by caricatures of pQCD calculations; see Eq. (1) of Ref. [13]. 
With Eq. (5.4) in hand, the form factors can be obtained in closed form. The results 
are 

2 

Fnn = 1 - e x y/nxErfc(x), (5.5) 
F 2)N = -x 2 + -e x2 ^xErfc(x) + e x " ^ x s Erfc(x) . (5.6) 
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where x = ^^j 1 and Erfc(x) is the complementary error function. We use Eq. (5.3) 
to calculate the cross section . Fig. 8 shows the finite sized result obtained using 
Eqs. (5.5) and (5.6), labelled by OBOfs, compared with the zero-size OBO\. We see 
that including finite size, via Eq. (5.4) has a very small effect. 

VI. N* PRODUCTION 
We can also apply our machinery to examine iV* production. Very little changes 
in the formalism if a PLC is formed and if the b 2 form of the wave-packet nucleon 
interaction is used. 

Let us first consider the DWBA for iV*'s. The N* experiences its own optical 
potential as it moves through the nucleus. We can obtain the result for the distorted 
wave by neglecting states other than the N* , e.g. make U diagonal in the quark space. 
The result is 

Vdwba(R) = e^ z e-^f-~ dZlPi§ ' Zl) = e ^ z e x ^ B ' z \ (6.1) 

where p2 the momentum of the outgoing excited state is close to the photon momentum 
q. Note the appearance of the factor 3 in the exponent of Eq. (6.1). This factor has a 
large effect, causing terms in which the N* decays to be more important than those in 
which the N* propagates. 

Using previous notation, we have calculated OBO\, OBO\, LEE{, LEE%, and 
EA*. The superscripts serve to remind us that we are calculating N* terms. 

We first show the OBO* results. Recall, V(R) = e ip2Z [$>° BO * (R) + V° BO * (R) + 
V$ BO *(R) + ---]. 

^ BO *(^) = l, (6.2) 

r z 

y° BO *(R) = --\ dZ lP (B, Z x )(3 - e l{pn ~ P2){z - Zl) - 2e lipi ~ P2){z - Zl) ) , (6.3) 
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*% BO \R) = (|) 2 f Z dZ 1 p(B,Z 1 ) r dZ 2 p(B,Z 2 )(9 

L J -oo J -oo 

_i_ 4 e i(P4-P2)(Z 1 -Z 2 ) _|_ e i(po-P2)(Z 1 -Z 2 ) _ e i(po-p 2 )(Z-Z 2 ) 

— 3e i ( po_P2 )( z-Zl ) — g e i (P4-p2)(^--Z'i) _ iQ e i(pi-P2)(z - z 2 ) 
+ Q e i( .Pe-Pi)( z - z ^+ i (p-i-P'2)(z-z 2 )y 

The OBO* curves are shown on Fig. 9 with the same notation as before. The lower 
solid line is the full distorted wave Born approximation, denoted by DWBA* ; the upper 
solid line is the cross section calculated using only the first order term of the DWBA and 
is labelled DWBAl. The difference between those two curves is much bigger than in 
the nucleon case (the relevant lines there would correspond to DWBA = .54, DWBAi 
= .48) This difference is simply due to the presence of the extra factor of 3 in Eq. (6.1). 
The resonance- like bump at lower energies is a false effect. Because the first order term 
starts out so big, which we see is because of the large value produced by the first order 
DWBA, the tendency of the phases in Eqs. (6.2-6.4) is to increase the cross section at the 
lower energies. It is not surprising that the order by order calculation fails here. Indeed, 
we learned from our experience with the nucleon that the OBO curves were unreliable 
at low energies. The fact that the exponent of the DWBA is large in magnitude only 
accentuates that failure. On these grounds, we expect the LEE* and the EA* curves, 
which have many more higher order effects, to do much better. 

Of course, the EA* is calculated by simply exponentiating the OBO*, Eq. (6.3). 
The LEE* results for the distorted waves are shown below. 

^ EE " (R) = e x ^ z \ (6.5) 
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-2e i{p4 - p2 ^ z - Zl) J, (6.6) 

Z Z / 

^ EE *(R) = (L) 2 e *3,z) f dZ 1 p(B,Z 1 ) f 1 dZ 2 p(B,Z 2 )U + 4e^- p ^- z ^ 

J — oo J — oo \ 

_|_ e i(Po-P2)(Z 1 -Z 2 ) _ 2 e i (P0-P2)(Z-Z 1 ) _ ^ e i(p 4 -p 2 )(Z-Z 1 ) 

_ 8e i(p4-P2)(Z-Z 2 ) _|_ Q e i( P 6-P4)(Z-Z 1 )+i( P4 - P2 )(Z-Z 2 ) 

Here, X(B,Z) is defined in Eq. (3.11). The LEE\, LEEo*, and EA* are shown on 
Fig. 10 as the dotted, dashed, and solid lines, respectively. As expected, there is good 
agreement among all three curves, especially below — Q 2 ~ 20 GeV 2 . 

It is curious that the LEE* works as well as it does since now the 0th order term 
does not reproduce the DWBA*. Therefore, one might think that instead of defining 
tt- = 1+A6 2 it would be more appropriate to define the operator Ab 2 by the assignment 
£r = 3 + Ab 2 and then use Eqs. (3.12 - 3.13). Calculations using this "bar" method 
show that the that convergence is better in the LEE* scheme than in the LEE*. This 
is due to the preference, noted above, that iV* decay is preferred to iV* propagation. 

Note that in this iV* work, we have taken the ejected wavepacket to be in a PLC. 
Qualitative changes in the predictions could occur if this is assumption is not correct [11]. 

We now summarize our results for the iV* case. 

1. The curves OBO\ and OBO*. do not really start to converge well at the energies 
shown. However, starting at about — Q 2 IGGeV 2 , the two curves are not far 
apart. 

2. The curves LEE^ and LEE\ show good agreement for energies less than — Q 2 « 
20GeV 2 . Above that, the LEE\ is not accurate. The fact that the LEE 2 agrees 
so well with the EA* starting from — Q 2 8GeV 2 suggests that the correct answer 
is given by these two curves above this energy. 

3. Below -Q 2 ~ 8GeV 2 the LEE^ agrees extremely well with the EA* . This indicates 
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that these two curves give the correct answer for this energy range. Also, this leads 
us to the conclusion that the shoulder which appears at low energies in the LEE\ 
at low energies is not real. 
4. Our experience in the nucleon sector taught us that the LEE and the EA were the 
best approximations. This conclusion remains true in the case of AT*'s too. 

VII. SUMMARY 

We have derived approximation schemes that are accurate at low (LEE, Sect. IIB) 
and high (HEE, Sect. IIC) energies. Their regimes of accuracy overlap, so that methods 
for calculations at any energy exist. The applications here use simplified interactions 
and baryon wavefunctions. The methods are more general and allow calculations using 
more realistic interactions and wavefunctions. 
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FIGURE CAPTIONS 

Fig. 1. Order-by-order approximation scheme. Ratios of cross sections for the 12 C(e 7 e'p) 
reaction. Dotted- OBOi. Dashed-OB0 2 . Solid- OB0 3 . 

Fig. 2. The low energy expansion. Dotted-LEEi. Dashed-LEE 2 . Solid-LEE 3 . 

Fig. 3. The high energy expansion. 

Fig. 4. The exponential approximation. 

Fig. 5. Comparing different approximation schemes. 

Fig. 6. Use of the interaction of Eq. (4.1), Ref. [15]. 

Fig. 7. The effect of a finite sized wave-packet at very high Q 2 . Size ~ j^. 

Fig. 8. Energy dependence of the finite-size effects, A ~ Q. 

Fig. 9. Order-by-order approximation for quasi-elastic N* production. 

Fig. 10. Low energy expansion for quasi-elastic N* production. 
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